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ABSTRACT 

The basic characteristics of the thermoactivated vortex mobility in Y 1 Ba 2 Cu 3 0 7 are 
determined by measurement of the kinetics of magnetization in two time regimes . 
The analysis of the kinetics of the approach of the equilibrium results in the 
activation energy, while the measurement of the log- creep rate allows 
determination of the activated moment. It is shown that the movement of vortices 
can be regarded as the diffusion process. 

INTRODUCTION 

The problem of vortex pinning is of great importance for the application of oxide 
high Tc superconductors. In comparison with pinning of dislocations, interface 
boundaries or domain walls, the pinning problem for vortices is much more complex 
because the density of vortices is large and adequate to produce strong 
interactions between the vortices in many practical important cases. 


In accordance with the effects of vortex-vortex interactions on pinning 
two different cases can be considered. The first one corresponds to strong 
pinning centers with the spacing much larger than the distance between vortices. 
In this case the system of vortices can be considered an equilibrium phase 
(crystal , glass or melt) , and the interaction of this phase with randomly 
distributed pinning obstacles can be analyzed. A classic example of such an 
approach is the Labusch theory of pinning for an elastic vortex lattice. 
Another example is Nelson's melt theory, 2 which pretends to the explanation of 
the very weak pinning at temperatures above the reversibility temperature. Note 
that before melting, a great number of dislocations must be formed in the vortex 
crystal, and the real problem of the dynamics of the vortex medium in this case 
is the problem of movement and pinning of dislocations in the vortex lattice. 
In other words the trend to melting has to be manifested by a transition from a 
laminar flow of the vortex medium to a turbulent one. This problem seems very 
important in the case of pinning by interfaces or particles with spacing larger 
than a distance between vortices. 3 

However, the problem of intrinsic pinning in oxide superconductors may be 
closer to reality in the other case, where the individual vortices are pinned by 
the obstacles, and the distance between the obstacles is comparable with the 
distance between vortices . Collective effects of the vortex-vortex interactions 
are important in this case also, but we can expect that they manifest themselves 
in the renormalization of parameter values rather than in the alteration of the 
physical picture of the vortex motion. This case of small vortex density is under 
consideration in the paper. Our goal is to obtain basic characteristics of vortex 
mobility from the study of the magnetization of zero-field cooled samples. 
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THEORETICAL CONSIDERATION 

The penetration of magnetic flow into a superconductor during its magnetization 
in an external magnetic field, Fig. 1, is similar to the flux flow diffusion , * > 5 
and is described by the generalized diffusion equation: 

db/dt - V(D eff VB) , ( 1 ) 

where p eff «=vB/VB is the effective flux diffusion coefficient, B is the local 
induction, VB is the gradient of the induction, and v is the flux velocity. 

The equation (1) is correct if there is no nucleation, multiplication 
or annihilation of vortices inside the sample, i.e. there are no weak links 
within a superconductor. The flux velocity for the thermoactivated movement of 
a vortex through an array of pinning obstacles can be determined as 


v - v c exp(-U/kT)sinh(aVB/kT) . (2) 

This equation is obtained under the assumption that the energy barrier AF, 
which has been overcome for a vortex displacement, is disturbed under the driving 
force (VB) , and it is an analytical function of VB at VB-O :*- 6 


AF ± (VB) - AF(0) T aVB, AF(0) - U, ( 3 ) 

where + corresponds to displacements with VB, - to the opposite direction, U 
is the activation energy, and a is the activated moment, a-mij, ra-^ 0 ni 2 , where 
m is the magnetic moment of the vortex or vortex bundle, ^ o -2,7xl0 ' 7 Gem 2 is the 
flux of one vortex, n is the number of the vortices in the bundle, i, is the 
displacement of the bundle or vortex at an activation event, i 2 is the effective 
length of a displaced section of the bundle or vortex, factually the distance 
between the nearest pinning points along the vortex, (Fig. 2). v 0 in Eq.(2) is 
the limiting velocity which equals i^, where v is the frequency of vibration of 
the vortex section, i 2 . Assuming i ■'~(£/i 2 )*' at . where £ is the coherence length or 
vortex "radius" (-10 cm), is the atomic vibration frequency - 10 13 s _1 , we get 
for v o “£i / at “10 cm/s, a value on the order of the sound velocity. 


If the driving force is larger than the critical force, VB>VB cr -U/a, the 
movement of the vortex is overbarrier. The vortex movement is not described by 
Eq. (2) because it is controlled not by the frequency of fluctuation depinning 
events but by the rate of the dissipation of the vortex energy. In this case 
v ”(l/ 7 )VB, where 7 is the viscosity . 7,8 


The vortex velocity, v(VB) , and the corresponding diffusion coefficient, 
D eff(? B) , are given in Fig. 3. There are three regions for which the solutioA 
or Eq . (1) can be obtained analytically. 

The first one corresponds to the distribution of magnetic induction with 
uR Cr ’ ThlS 1S thS S ° called flux flow re g io n when D eff =D f -B / 7 is independent 


The second region corresponds to the near critical state, when aVB/kT»l 
and sinh(aVB/kT)=exp(aVB/kT) . In this region the total magnetization of ’a 
superconductor, or its average induction, is proportional to the logarithm of 
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time : A,e 

M(t) - P lnt, P “ dM/dlnt « kTR/(3AF/3VB) - (kTR/a(R)), 

where R is the radius of the sample, a(R) is the activated moment at its surface 
where the induction B(R) is closed to the equilibrium one, B 0 . 


The third region corresponds to the constant D eff =D at the small VB 

again: 


D - ( B 0 ctv 0 /kT ) exp ( - U/kT ) - D 0 exp( -U/kT) , D 0 - B 0 av 0 /kT 


(5) 


This regime of ordinary diffusion flow was named the thermally assisted flux flow 
(TAFF) . 8 


It is important to note, that all of these regimes can be observed by 
measuring the magnetization as a function of time of the zero- field pooled 
superconductor sample. At the first moment after applying the magnetic field, VB 
is great and is larger than the critical value. Unfortunately, this stage is 
difficult to study because of nonstationary effects while the field is turned 
on. The sample then goes through the near critical state, where linear behaviour 
of magnetization vs. logarithm of time has been expected. At last, in the final 
stage, the magnetization approaches to the equilibrium saturation, Mo(° r o) > a 
t>t* , * and the kinetics of the magnetization is described by the equation: 


M(t) - ^ + M ' exp ( - t/r ) , 


( 6 ) 


where M' is a constant. The relaxation time, r , is connected with the diffusion 
coefficient, D(B 0 ) , by the relation r - R 2 /D, where R is the sample radius. It 
is necessary to stress that this near equilibrium asymptotic behaviour does not 
depend on the distribution of the local induction at the moment t*, when the 
exponential stage begins , provided the induction is closed to B 0 (RVB/B 0 « ) . 

Thus the study of the kinetics of the magnetization at the different stages 
gives the opportunity to find the basic characteristics of the vortex mobility, 
a and U, without any additional assumption about the distribution of the 
induction within the sample. The results of such a kind of study are represented 
below. 


SAMPLES AND EXPERIMENTAL METHOD 


Three different Y^azCuaO; samples with sharp transitions at Tc-93K were under 
investigation: a single crystal, a sintered polycrystal and a sample of 

quench, melt and growth (QMG) material. 10 The QMG sample is a bar of 0. 9x0. 9x5 .^8 
mm, which was cut from a grain of a polycrystal. The sample has no weak links, 
making it possible to determine its effective radius, R (R-0.45 mm) . Our single 
crystal sample contains some fraction of normal phase. So, the single crystal as 
well as sintered samples were not good for quantative estimates. 


The samples were cooled in zero magnetic field. The magnetic field was then 
applied, and the magnetization was measured as a function of time using either 
a vibrating sample magnetometer, (VSM) , or a SQUID system. The magnetic field was 
directed along c axes of the single crystal sample and perpendicularly to c axes 
of the QMG sample. 
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EXPERIMENTAL RESULTS 


Approach to Equilibrium Magnet ization , M ,. Diffusion Coefficient D Magnetization 
versus time, M(t), curves were measured at rather high temperatures (70 and 80K) 
lon S exposures (up to 340 ksec) at relatively weak magnetic fields (2 and 
4 kOe) so the density of vortices was small up to approaching the equilibrium 
magnetization. M(t) data for all the samples at 70K and 2 kOe are plotted in 
curves (a) of Figures 4,5 and 6 for the QMG, single crystal and sintered 
samples respectively. M vs. log- time data are plotted in curves (b) . Curves 
obtained for the other temperature and field setting have the same shape as in 
these Figures: for all the curves there are log- time dependences M(t), and at 
much longer times, t>t,* there are stages of the approach to the equilibrium 

behaviour tl0n> M °* ^ these stages the magnetization shows the exponential 


ln[M(t)/M 0 - 1] - ln(M'/Mo) “ -t /r , (7) 

where ln(M'/M 0 ) are constants. In curves (c) of Fig. 4, 5 and 6, where the 
calculated data of ln[M(t)/M 0 -l] vs. time for all the samples are plotted, it can 
be seen that at t>t* there are linear dependences of ln[M(t)/M 0 -l] on time. The 
result of the determination of r depends strongly on the choosen value of M and 
asymptotic method was used 12 for the correction the values of m!' For 
the QMG sample we obtained the values of M.-3.87 (emu/gr) at 70K and 2 kOe, and 
M o~- 2,5 (emu/gr) at 70K and 4 kOe. 

The values of the relaxation time, r, were obtained as a slope of the 
linear portion of the curves (c) , and the diffusion coefficients, D-R 2 /r, where 

?n£ S the Kt ffeC a 1 n e / a -, dl ;n S 9 0f p he Saraple ’ were calculated. For the QMG sample at 
70K we obtained D-6 . 7x10 cm 2 /s for 2 kOe, and D=7.8xlCT 9 cm 2 /s for 4 kOe. 


Activated Moment , a, Activated Length, l. and Intervortex Snacinr * The 
logarithmic creep rate P-dM/dlnt for the QMG and sintered samples was measured 
in the temperature interval 6-80K. 10 The results are given in Fig. 7, a and in 
g. 7,b, respectively. For the QMG sample we calculated the values of the 
activated moment, a (using the relation o-kTR/P, Eq. (4)). Derived values of 
a/4> 0 , where <f> 0 is the flux of one vortex, for the QMG sample are plotted in Fig. 


Values of the average activated length, i-W^) 1 ' 2 , are displayed in Fig. 
9 together with values of the distance between vortices, a-(B /6 ) 1/2 . B was 
obtained from measuring the equilibrium magnetization, M 0 : Bq-^ttMo+H. ° 

Activation Energy, U, and Pre- exponential Factor, P 3 . Obtaining r and a from the 
p^ 1CS u° magnetization, it is possible to estimate the activation energy 

of the thermoactivated movement of the vortex bundle, U, from the Eq. (5): 


U/kT = ln(D/D 0 ) = ln(rv 0 B 0 a/kTR 2 ) , 


( 8 ) 


Assuming that v o -10 5 cm/s, we obtained, with log-accuracy by v U/kT=33 U^O 19 
eV, and the pre-exponential factor D o -B o av o /kT~6xl0 5 cm 2 /s for °the 
70K and 2 kOe . For 70K and 4 kOe U-0.2 eV, D o =*4xl0 6 cm 2 /s. 


QMG sample at 
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DISCUSSION 


The exponential approach to equilibrium is evidence of the linear dependence 
between the vortex velocity and the driving force and supports the assumption of 
Eq .(3). Taking into account the relation between D G ff and resistance, P * P ™ 
D/c 2 , where c is the velocity of light, we have to admit that at 70 and 80K and 
H-2 and 4 kOe the QMG sample of YBa 2 Cu 3 0 7 manifests small but finite 
resistivity. From the point of view of vortex mobility in this stage, the vortex 
medium is rather "melt” than "glass" (in terms of M. Fisher's theory 13 ) or 
"crystal" with collective pinning. 14 For both these models AF must be the 
divergent function of VB at V-0: AF-l/VB* 3 , £>0 for the vortex glass state, or 
-InVB for the collective pinning model. 

By "vortex melt" , we imply that the effects of correlations between the 
individual thermoactivated depinning events are small. 15 To prove this suggestion 
it is necessary to consider the pre- exponential factor, 0 o , in Eq. (5) in more 
detail, rewriting it through the activated length, i, and the intervortex 
distance, a , as 

D 0 = (^ 0 2 /kT)(i 2 /a 2 )v 0 =* Lv 0 (j e 2 /* 2 ). < 9 ) 

The term, ^ 0 2 /kT»L, has the dimension of length (e.g. l/=*4 cm at 70K) . From Fig. 
10 it is clear that up to 70K at 2 kOe and 60K at 4 kOe, i<a, i.e. a bundle 
consists of only one vortex, and £ is the weak function of temperature. Then £ 
rises rapidly, i.e. the number of vortices in the bundle increases. For n-1, 
X-(i 1 i 2 ) 1/2 , that is, £ is close to an average distance between pinning points. 

At 7 OK and higher temperatures, all the pinning centers are occupied by 
vortices, and the pinning -depinning movements of vortices with the driving force 
VB can be interpreted as back-directed diffusion of pinning points in vortex 
medium. Indeed, the velocity of the point flow is determined by the Einstein 
diffusion relation: 


v p - (D p /kT) V/i = (Dp/kT) (aVB/a) , 


( 10 ) 


where D p -a z yexp( -U/kT) is the point diffusion coefficent, and V/j is the gradient 
of chemical potential of the pinning points in the magnetic field. The point 
velocity, v p is equal numerically to the vortex velocity, v: 

-v p _ v - DVB/B (11) 

(the relation between D and v arises from the definition of D eff in Eq. (1)) . If 
we take into account that av— v 0 , it follows from Eqs. (10) and (11) that the 
vortex diffusion coefficient, D, is determined by Eq, (5). At i=*a, 
D o «Lv o “10 6 cm 2 /s: that is in accordance with the experimental determination of D 0 
through measuring a and B 0 . 


It is clear that the penetration of the flux flow into a sample during 
magnetization can be considered in terms of the diffusion theory for polycrystal 
samples with weak grain boundaries too. For this case it would be necessary to 
bring into consideration the vortex diffusion along the grain boundaries as well 
as the bulk grain diffusion. Of course the diffusion description may be not 
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complete if the nucleation of vortices is possible inside "bad” samples. 
CONCLUSION 

1. It is shown that approach to equilibrium magnetization with time of the zero- 
field cooled QMG, single crystal and sintered samples of YBa 2 Cu 3 0 7 proceeds 
according to the exponential law. The equilibrium values of the induction, 
relaxation time and vortex diffusion coefficient are determined for the QMG 
sample at 70 and 80K and 2 and 4 kOe . 

2. Log-creep rate was measured at 2 and 4 kOe for a set of temperatures between 
6 and 80K. From these data the activated moment a was found. The analysis of a 
shows that an activated vortex bundle contains only one vortex at the 
temperatures and fields under consideration. An activation length, i - (a/^ Q ) 1/2 , 
is almost constant up to 70 or 80 K. Beyond these temperatures the distance 
between vortices becomes smaller than i. Up to these temperatures £ is likely the 
average distance between effective pinning points. At temperatures higher than 
7 OK at field 2 kOe and temperatures -75K at 4 kOe the number of vortices in a 
bundle grows, as does £. 

3. Using B 0 and D, which are obtained from the approach to the equilibrium , and 
a from the log-creep stage, the activation energy, U, has been determined for the 
QMG sample. 

4. The character of vortex movement during magnetization with time at 70 and 80K 
and 2 and 4 kOe is In good agreement with the diffusion picture, which does not 
take into account the effects due to the correlation between individual 
thermo activation events . 
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(a) (b) 



Fig. 1. Kinetics of magnetization of a zero-field cooled sample. 

a. Equilibrium magnetization dependence on magnetic field. 

b. Time dependence of magnetization. 

c. Distribution of magnetic induction. 

0 - the initial diamagnetic state, 1-2- the near critical state (VB<VB cr ) , 
3 - the near equilibrium state, 4 - the equilibrium state (VB-0, B-B 0 ) . 




Fig. 2. Thermoactivated movement of a vortex, a. A vortex section, i 2 , pinned 
by obstacles, b. Energy as function of a vortex displacement. 

Fig. 3. Vortex velocity, v, and vortex diffusion coefficient, D eff . 
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M («rau/gr) ^ ^ M (eBU^gr) ^ ln(M/Mo.l) * M (eau/gr) * M (eou/gr) 



In (tl»e, sec) 


In (ti»« t sec) 













dM/d(ln t) [emu/g] p = d M/d(ln t) [emv/g] 



Fig. 7. F 


Fig. 7. Logarithmic rate of magnetization, P. 

a. the QMG sample, b. the sintered sample. 

Fig. 8. Temperature dependence of a/<f> 0 for the QMG sample at 
Fig. 9. Temperature dependence of average activated length, 
spacing, a, for the QMG sample at 2 and 4 kOe. 






